Reduced genetic diversity and the success of the invasive peacock bass (Cichliformes: Cichlidae) / Diversidade genética reduzida e o sucesso da espécie invasora tucunaré (Cichliformes: Cichlidae)

Abstract Several species of Cichla successfully colonized lakes and reservoirs of Brazil, since the 1960's, causing serious damage to local wildlife. In this study, 135 peacock bass were collected in a reservoir complex in order to identify if they represented a single dominant species or multiple ones, as several Cichla species have been reported in the basin. Specimens were identified by color pattern, morphometric and meristic data, and using mitochondrial markers COI, 16S rDNA and Control Region (CR). Overlapping morphological data and similar coloration patterns prevented their identification using the taxonomic keys to species identification available in the literature. However, Bayesian and maximum likelihood from sequencing data demonstrated the occurrence of a single species, Cichla kelberi. A single haplotype was observed for the 16S and CR, while three were detected for COI, with a dominant haplotype present in 98.5% of the samples. The extreme low diversity of the transplanted C. kelberi evidenced a limited number of founding maternal lineages. The success of this colonization seems to rely mainly on abiotic factors, such as increased water transparency of lentic environments that favor visual predators that along with the absence of predators, have made C. kelberi a successful invader of these reservoirs.

Resumo Muitas espécies de Cichla colonizaram com sucesso lagos e reservatórios do Brasil desde os anos 1960, causando graves prejuízos à vida selvagem nesses locais. Neste estudo, 135 tucunarés foram coletados em um complexo de reservatórios a fim de identificar se representavam uma espécie dominante ou múltiplas espécies, uma vez que diversas espécies de Cichla foram registradas na bacia. Os espécimes foram identificados com base na coloração, dados morfométricos e merísticos, e por marcadores mitocondriais COI, 16S rDNA e Região Controle (RC). A sobreposição dos dados morfométricos e o padrão similar de coloração impediram a identificação utilizando as chaves de identificação disponíveis na literatura. Entretanto, as análises bayesiana e de máxima verossimilhança de dados moleculares demonstraram a ocorrência de uma única espécie, Cichla kelberi. Um único haplótipo foi observado para o 16S e RC, enquanto três foram detectados para o COI, com um haplótipo dominante presente em 98,5% das amostras. A baixa diversidade nos exemplares introduzidos de C. kelberi evidenciou um número limitado de linhagens maternas fundadoras. O sucesso da invasão parece depender de fatores abióticos, como a maior transparência da água de ambientes lênticos que favorece predadores visuais que, atrelado à ausência de predadores, fez do C. kelberi um invasor bem-sucedido nesses reservatórios.

Isolation and characterization of bacteria associated with silkworm gut under antibiotic-treated larval feeding / Isolamento e caracterização de bactérias associadas ao trato intestinal do bicho-da-seda sob alimentação larval tratada com antibióticos

Abstract The impact of antibiotics on growth, cocoon production was assessed in addition to isolation and characterization of bacteria associated with silkworm gut of infected larvae. Larval rearing was maintained at recommended conditions of temperature and humidity. Silkworm larvae showing abnormal symptoms were collected from the control group and dissected for gut collection. Bacteria were isolated from the gut content by spreading on agar plates and incubated at 37 °C for 48 hrs. Bacterial identification and phylogenetic analysis were carried out by 16S rRNA gene sequencing. The isolated bacteria were subjected to antimicrobial susceptibility test (disc diffusion methods) by using Penicillin (10 µg/mL), Tetracycline (30 µg/mL), Amoxicillin (25 µg/mL), Ampicillin (10 µg/mL), and Erythromycin (15 µg/mL). All isolated strains showed positive results for the catalase test. We isolated and identified bacterial strains (n = 06) from the gut of healthy and diseased silkworm larvae. Based on the 16S rRNA gene sequence, isolated bacteria showed close relation with Serratia, Bacillus, and Pseudomonas spp. Notably, 83.3% of strains were resistant to Penicillin, Tetracycline, Amoxicillin, Ampicillin, and Erythromycin but 16.6% showed antibiotic susceptibility to the above-mentioned commonly used antibiotics. Silkworm larvae fed on penicillin-treated leaves showed significant improvement in larval weight, larval length, and cocoon production. Significantly higher larval weight (6.88g), larval length (5.84cm), and cocoon weight (1.33g) were recorded for larvae fed on leaves treated with penicillin as compared to other antibiotics. Isolated bacterial strains showed close relation with Serratia spp., Bacillus spp. and Pseudomonas spp.

Resumo O impacto dos antibióticos no crescimento e na produção do casulo foi avaliado, além do isolamento e caracterização das bactérias associadas ao intestino de larvas infectadas do bicho-da-seda. A criação das larvas foi mantida nas condições recomendadas de temperatura e umidade. As larvas do bicho-da-seda com sintomas anormais foram coletadas do grupo controle e dissecadas para coleta do intestino. As bactérias foram isoladas do conteúdo intestinal por espalhamento em placas de ágar e incubadas a 37° C durante 48 horas. A identificação bacteriana e a análise filogenética foram realizadas pelo sequenciamento do gene 16S rRNA. As bactérias isoladas foram submetidas a teste de sensibilidade antimicrobiana (métodos de difusão em disco) com penicilina (10 µg / mL), tetraciclina (30 µg / mL), amoxicilina (25 µg / mL), ampicilina (10 µg / mL) e eritromicina (15 µg / mL). Todas as cepas isoladas apresentaram resultados positivos para o teste da catalase. Isolamos e identificamos cepas bacterianas (n = 06) do intestino de larvas de bicho-da-seda saudáveis e doentes. Com base na sequência do gene 16S rRNA, as bactérias isoladas mostraram estreita relação com Serratia, Bacillus e Pseudomonas spp. Notavelmente, 83,3% das cepas eram resistentes a penicilina, tetraciclina, amoxicilina, ampicilina e eritromicina, mas 16,6% mostraram suscetibilidade aos antibióticos comumente usados mencionados acima. As larvas do bicho-da-seda alimentadas com folhas tratadas com penicilina apresentaram melhora significativa no peso larval, comprimento larval e produção de casulo. Peso larval significativamente maior (6,88g), comprimento larval (5,84cm) e peso do casulo (1,33g) foram registrados para larvas alimentadas com folhas tratadas com penicilina, em comparação com outros antibióticos. Cepas bacterianas isoladas mostraram estreita relação com Serratia spp., Bacillus spp. e Pseudomonas spp.

Genetic variations and phylogenetic relationship of genus Uromastyx from Punjab Pakistan / Variações genéticas e relação filgenética do gênero Uromastyx de Punjab, Paquistão

Abstract During the present study, specimens were collected from selected sites of Cholistan desert and Kalabagh Game Reserve, Punjab province, Pakistan. Each captured specimen was tagged with voucher number and morphometric measurements were taken. The average snout to vent length was 172.559±1.40 mm and average weight was 92.1±1.30 g. The DNA of Uromastyx hardwickii was amplified and sequenced using 16S rRNA primer set. The obtained DNA sequence has shown reliable and clear species identification. After trimming ambiguous bases, the obtained 16S rRNA fragment was 520 bp while 16S rRNA fragments aligned with closely matched sequence from NCBI comprised of 510 bp. Closely matched sequences of genus Uromastyx were retrieved from NCBI in blast searches. Neighbour-joining tree of genus Uromastyx was constructed based on p-distance using MEGA X. The mean intraspecific variation was 0.095±0.01 while intraspecific variation was ranging from 0-1%. Similarly, interspecific variation of Uromastyx hardwikii with Saara asmussi, Uromastyx alfredschmidti, Uromastyx geyri, Uromastyx thomasi, Uromastyx alfredschmidti was 0-12%, 0-19%, 0-19%, 0-20%, 12-19% respectively. The newly produced DNA was submitted to NCBI and accession number was obtained (MW052563.1). Results of current study provided information about the molecular and morphological identification of Genus Uromastyx. In our recommendation, comprehensive molecular based identification of Pakistan's reptiles is required to report any new or subspecies from country.

Resumo Durante o presente estudo, os espécimes foram coletados em locais selecionados do deserto do Cholistan e da Reserva de Caça de Kalabagh, província de Punjab, Paquistão. Cada espécime capturado foi etiquetado com o número do comprovante e medidas morfométricas foram realizadas. O comprimento médio do focinho à cloaca foi de 172,559 ± 1,40 mm, e o peso médio foi de 92,1 ± 1,30 g. O DNA de Uromastyx hardwickii foi amplificado e sequenciado usando o conjunto de primer 16S rRNA. A sequência de DNA obtida mostrou identificação de espécies confiável e clara. Após o corte de bases ambíguas, o fragmento de rRNA 16S obtido tinha 520 pb, enquanto os fragmentos de rRNA 16S alinhados com a sequência próxima do NCBI composta por 510 pb. Sequências semelhantes do gênero Uromastyx foram recuperadas do NCBI em pesquisas de explosão. A árvore de união de vizinhos do gênero Uromastyx foi construída com base na distância-p usando MEGA X. A variação intraespecífica média foi de 0,095 ± 0,01, enquanto a variação intraespecífica foi de 0-1%. Da mesma forma, a variação interespecífica de Uromastyx hardwikii com Saara asmussi, Uromastyx alfredschmidti, Uromastyx geyri, Uromastyx thomasi, Uromastyx alfredschmidti foi de 0-12%, 0-19%, 0-19%, 0-20%, 12-19%, respectivamente. O DNA recém-produzido foi submetido ao NCBI e o número de acesso foi obtido (MW052563.1). Os resultados do estudo atual forneceram informações sobre a identificação molecular e morfológica do Gênero Uromastyx. Em nossa recomendação, a identificação de base molecular abrangente de répteis do Paquistão é necessária para relatar qualquer nova ou subespécie do país.

Diversity of culturable endophytic fungi vary through time in Manihot esculenta Crantz / A diversidade de fungos endofíticos cultiváveis ​​varia ao longo do tempo em Manihot esculenta Crantz

Abstract Endophytic fungi are a ubiquituos group that colonize all plant species on earth. Studies comparing the location of endophytic fungi within the leaves and the sampling time in Manihot esculenta Crantz (cassava) are limited. In this study, mature leaves of M. esculenta from Panama were collected in order to compare the cultivable diversity of endophytic fungi and to determine their distribution within the leaves. A total of one hundred sixty endophytes belonging to 97 species representing 13 genera and 8 morphospecies determined as mycelia sterilia that containing 63 isolates were isolated. Cladosporium, Nigrospora, Periconia, and mycelia sterilia 1 and 3 were the most predominant isolated endophytes. We detected that endophytes varied across the sampling time, but not amongst locations within leaves. The endophytes composition across sampling and the location of endophytes within leaf was similar, except for Periconia and mycelia sterilia 3 and 7. The data generated in this study contribute to the knowledge on the biodiversity of endophytic fungi in Panama, and establish the bases for future research focused on understanding the function of endophytes in M. esculenta crops.

Resumo Os fungos endofíticos são um grupo ubiquituo que colonizam todas as espécies de plantas na terra. Os estudos que comparam a localização dos fungos endofíticos dentro das folhas de Manihot esculenta Crantz (mandioca) e o tempo de amostragem são muito escassos. Neste estudo, folhas maduras de M. esculenta foram coletadas do Panamá com a finalidade de comparar a diversidade cultivável de endófitos e determinar sua distribuição dentro das folhas. Um total de 170 endófitos foram isolados de 97 espécies que representam 13 gêneros e 8 morfoespécies determinadas como micélios esterilizados contendo 63 isolados. Os fungos Cladosporium, Nigrospora, Periconia e mycelia sterilia 1 e 3 foram os isolados mais predominantes. Também detectamos que os endófitos variaram ao longo do tempo de amostragem, mas não entre os locais dentro das folhas. A composição de endófitos na amostragem e localização de endófitos dentro da folha foi semelhante, exceto para Periconia e mycelia sterilia 3 e 7. Os dados gerados neste estudo contribuem para o conhecimento da biodiversidade de fungos endofíticos no Panamá e estabelecem as bases para pesquisas sobre o entendimento da função de endófitos em culturas de M. esculenta.

Genomic constellations of RVA detected in Brazil from 1986 to 2016: a temporal and geographical distribution and occurrence of reassortments / Constelações genômicas de RVA detectadas no Brasil de 1986 a 2016: uma distribuição temporal e geográfica e ocorrência de rearranjos

INTRODUCTION: Species A rotavirus (RVA) infections are a major cause of severe gastroenteritis in children of <5 years worldwide. In Brazil, before vaccination, RVA was associated with 3.5 million episodes of acute diarrheal disease per year. Due to the segmented nature of their genomes, rotaviruses can exchange genes during co-infections, and generate new virus strains and new reinfections. OBJECTIVE: To evaluate the genomic diversity of RVA isolated in Brazil in 30 years, between 1986 to 2016, to investigate possible changes in the frequency of genotype constellations before and after the implementation of the vaccine. METHODS: In total, 4,474 nucleotide sequences were obtained from the Virus Variation Database. Genomic constellation was compared, and the proportion of rotavirus genotypes was analyzed by time and geographic region. RESULTS: Our results showed that major known genotypes were circulating in the country during the period under analysis, with a prevalence of the G1P[8] Wa-like genotype, decreasing only in the period immediately after the introduction of the vaccine. Regarding the geographical distribution, most of our constellations, 62 (39.2%), and 50 (31.6%) were concentrated in the North and Northeast regions. Our analysis also showed the circulation of multiple strains during the periods when the DS-1-like and AU-1-like genotypes were co-circulating with the Wa-like genotype. CONCLUSION: Therefore, it is likely that inter-genogroup reassortments are still occurring in Brazil and so it is important to establish an efficient surveillance system to follow the emergence of novel reassorted strains that might not be targeted by the vaccine.

INTRODUÇÃO: As infecções por rotavírus A (RVA) são uma das principais causas de gastroenterite grave em crianças <5 anos em todo o mundo. No Brasil, antes da vacinação, o RVA estava associado a 3,5 milhões de episódios de diarreia aguda por ano. Devido à natureza segmentada de seus genomas, os rotavírus podem trocar genes durante as coinfecções, gerar novas cepas de vírus e novas reinfecções. OBJETIVO: Avaliar a diversidade genômica de RVA isolados no Brasil entre 1986 a 2016 para investigar possíveis alterações na frequência das constelações de genótipos antes e após a implantação da vacina. MÉTODOS: No total, 4.474 sequências de nucleotídeos foram obtidas do Banco de Dados de Variação de Vírus. A constelação genômica foi comparada e a proporção dos genótipos de rotavírus foi analisada por tempo e região geográfica. RESULTADOS: Nossos resultados mostraram que os principais genótipos conhecidos circulavam no país no período em análise, com prevalência do genótipo G1P[8] Wa-like, diminuindo apenas no período imediatamente após a introdução da vacina. Em relação à distribuição geográfica, a maioria das nossas constelações, 62 (39,2%) e 50 (31,6%), concentrava-se nas regiões Norte e Nordeste. Nossa análise também mostrou a circulação de cepas múltiplas durante os períodos em que os genótipos DS-1-like e AU-1-like estavam co-circulando com o genótipo Wa-like. CONCLUSÃO: Portanto, é provável que rearranjos inter-genogrupos ainda estejam ocorrendo no Brasil e por isso é importante estabelecer um sistema de vigilância eficiente para acompanhar o surgimento de novas cepas rearranjadas que podem não ser protegidas pela vacina.

Genome-wide analysis of AP2/ERF superfamily in Isatis indigotica / 中西医结合学报

OBJECTIVE@#AP2/ERF (APETALA2/ethylene-responsive factor) superfamily is one of the largest gene families in plants and has been reported to participate in various biological processes, such as the regulation of biosynthesis of active lignan. However, few studies have investigated the genome-wide role of the AP2/ERF superfamily in Isatis indigotica. This study establishes a complete picture of the AP2/ERF superfamily in I. indigotica and contributes valuable information for further functional characterization of IiAP2/ERF genes and supports further metabolic engineering.@*METHODS@#To identify the IiAP2/ERF superfamily genes, the AP2/ERF sequences from Arabidopsis thaliana and Brassica rapa were used as query sequences in the basic local alignment search tool. Bioinformatic analyses were conducted to investigate the protein structure, motif composition, chromosome location, phylogenetic relationship, and interaction network of the IiAP2/ERF superfamily genes. The accuracy of omics data was verified by quantitative polymerase chain reaction and heatmap analyses.@*RESULTS@#One hundred and twenty-six putative IiAP2/ERF genes in total were identified from the I. indigotica genome database in this study. By sequence alignment and phylogenetic analysis, the IiAP2/ERF genes were classified into 5 groups including AP2, ERF, DREB (dehydration-responsive element-binding factor), Soloist and RAV (related to abscisic acid insensitive 3/viviparous 1) subfamilies. Among which, 122 members were unevenly distributed across seven chromosomes. Sequence alignment showed that I. indigotica and A. thaliana had 30 pairs of orthologous genes, and we constructed their interaction network. The comprehensive analysis of gene expression pattern in different tissues suggested that these genes may play a significant role in organ growth and development of I. indigotica. Members that may regulate lignan biosynthesis in roots were also preliminarily identified. Ribonucleic acid sequencing analysis revealed that the expression of 76 IiAP2/ERF genes were up- or down-regulated under salt or drought treatment, among which, 33 IiAP2/ERF genes were regulated by both stresses.@*CONCLUSION@#This study undertook a genome-wide characterization of the AP2/ERF superfamily in I. indigotica, providing valuable information for further functional characterization of IiAP2/ERF genes and discovery of genetic targets for metabolic engineering.

Mitochondrial genome sequence characteristics and phylogenetic analysis of Schizothorax argentatus / 生物工程学报

Schizothorax argentatus that only distributes in the Ili River basin in Xinjiang is one of the rare and endangered species of schizothorax in China, thus has high scientific and economic values. In this study, the complete mitochondrial genome sequence of S. argenteus with a length of 16 580 bp was obtained by high-throughput sequencing. The gene compositions and arrangement were similar to those of typical vertebrates. It contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding region (D-loop). The nucleotide compositions were A (30.25%), G (17.28%), C (27.20%), and T (25.27%), respectively, showing obvious AT bias and anti-G bias. Among the tRNA genes, only tRNA-Ser(GCU) could not form a typical cloverleaf structure due to the lack of dihydrouracil arm. The AT-skew and GC-skew values of the ND6 gene were fluctuating the most, suggesting that the gene may experience different selection and mutation pressures from other genes. The mitochondrial control region of S. argenteus contained three different domains, i.e., termination sequence region (ETAS), central conserved region (CSB-F, CSB-E, CSB-D, and CSB-B), and conserved sequence region (CSB1, CSB2, and CSB3). The conserved sequence fragment TT (AT) nGTG, which was ubiquitous in Cypriniformes, was identified at about 50 bp downstream CSB3. Phylogenetic relationships based on the complete mitochondrial genome sequence of 28 Schizothorax species showed that S. argenteus had differentiated earlier and had a distant relationship with other species, which may be closely related to the geographical location and the hydrological environment where it lives.

Analysis of the chloroplast genome of Incarvillea younghusbandii Sprague / 生物工程学报

Incarvillea younghusbandii Sprague is a traditional tonic herb. The roots are used as herbal medicine for nourishing and strengthening, as well as treating postpartum milk deficiency and weakness. In this study, the chloroplast genome of I. younghusbandii was sequenced and assembled by the high-throughput sequencing technology. The sequence characteristics, sequence repeats, codon usage bias, phylogenetic relationships and estimated divergence time of I. younghusbandii were analyzed. The 159 323 bp sequence contained a large single copy (80 197 bp), a small single copy (9 030 bp) and two inverted repeat sequences (35 048 bp). It contained 120 genes, including 77 protein coding genes, 8 ribosomal RNA genes and 35 transfer RNA genes. AAA was the most frequent codon in the chloroplast coding sequence of I. younghusbandii. A total of 42 simple sequence repeats were identified in the chloroplast genome. Phylogenetic analysis revealed I. younghusbandii was mostly like its taxonomically close relative Incarvillea compacta. The divergence between I. younghusbandii and I. compacta was dated to 4.66 million years ago. This study was significant for the scientific conservation and development of resources related to I. compacta. It also provides a basic genetic resource for the subsequent species identification of the genus Incarvillea, and the population genetic diversity study of Bignoniaceae.

Characteristics and phylogenetic analysis of chloroplast genome of a new type of fruit Rubus rosaefolius / 生物工程学报

The genomic DNA of Rubus rosaefolius was extracted and sequenced by Illumina NovaSeq platform to obtain the complete chloroplast genome sequence, and the sequence characteristics and phylogenetic analysis of chloroplast genes were carried out. The results showed that the complete chloroplast genome of the R. rosaefolius was 155 650 bp in length and had a typical tetrad structure, including two reverse repeats (25 748 bp each), a large copy region (85 443 bp) and a small copy region (18 711 bp). A total of 131 genes were identified in the whole genome of R. rosaefolius chloroplast, including 86 protein coding genes, 37 tRNA genes and 8 rRNA genes. The GC content of the whole genome was 36.9%. The genome of R. rosaefolius chloroplast contains 47 scattered repeats and 72 simple sequence repeating (SSR) loci. The codon preference is leucine codon, and the codon at the end of A/U is preferred. Phylogenetic analysis showed that R. rosaefolius had the closest relationship with R. taiwanicola, followed by R. rubraangustifolius and R. glandulosopunctatus. The chloroplast genome characteristics and phylogenetic analysis of R. rosaefolius provide a theoretical basis for its genetic diversity research and chloroplast development and utilization.

Characteristics of the chloroplast genome of Dracaena marginata and phylogenetic analysis / 生物工程学报

Dracaena marginata is a widely cultivated horticultural plant in the world, which has high ornamental and medicinal value. In this study, the whole genome of leaves from D. marginata was sequenced by Illumina HiSeq 4000 platform. The chloroplast genome were assembled for functional annotation, sequence characteristics and phylogenetic analysis. The results showed that the chloroplast genome of D. marginata composed of four regions with a size of 154 926 bp, which was the smallest chloroplast genome reported for Dracaena species to date. A total of 132 genes were identified, including 86 coding genes, 38 tRNA genes and 8 rRNA genes. Codon bias analysis found that the codon usage bias was weak and there was a bias for using A/U base endings. 46 simple sequence repeat and 54 repeats loci were detected in the chloroplast genome, with the maximum detection rate in the large single copy region and inverted repeat region, respectively. The inverted repeats boundaries of D. marginata and Dracaena were highly conserved, whereas gene location differences occurred. Phylogenetic analysis revealed that D. serrulata and D. cinnabari form a monophyletic clade, which was the closest relationship and conformed to the morphological classification characteristics. The analysis of the chloroplast genome of D. marginata provides important data basis for species identification, genetic diversity and chloroplast genome engineering of Dracaena.

Chloroplast genomic characterization and phylogenetic analysis of Pellionia scabra / 生物工程学报

Pellionia scabra belongs to the genus Pellionia in the family of Urticaceae, and is a high-quality wild vegetables with high nutritional value. In this study, high-throughput techniques were used to sequence, assemble and annotate the chloroplast genome. We also analyzed its structure, and construct the phylogenetic trees from the P. scabra to further study the chloroplast genome characteristics. The results showed that the chloroplast genome size was 153 220 bp, and the GC content was 36.4%, which belonged to the typical tetrad structure in P. scabra. The chloroplast genome encodes 130 genes, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes in P. scabra. Among them, 15 genes contained 1 intron, 2 genes contained 2 introns, and rps12 had trans-splicing, respectively. In P. scabra, chloroplast genomes could be divided into four categories, including 43 photosynthesis, 64 self-replication, other 7 coding proteins, and 4 unknown functions. A total of 51 073 codons were detected in the chloroplast genome, among which the codon encoding leucine (Leu) accounted for the largest proportion, and the codon preferred to use A and U bases. There were 72 simple sequence repeats (SSRs) in the chloroplast genome of P. scabra, containing 58 single nucleotides, 12 dinucleotides, 1 trinucleotide, and 1 tetranucleotide. The ycf1 gene expansion was present at the IRb/SSC boundary. The phylogenetic trees showed that P. scabra (OL800583) was most closely related to Elatostema stewardii (MZ292972), Elatostema dissectum (MK227819) and Elatostema laevissimum var. laevissimum (MN189961). Taken together, our results provide worthwhile information for understanding the identification, genetic evolution, and genomics research of P. scabra species.

Identification and expression analysis of MADS-box gene family in Docynia delavayi (Franch.) Schneid / 生物工程学报

MADS-box gene family is a significant transcription factor family that plays a crucial role in regulating plant growth, development, signal transduction, and other processes. In order to study the characteristics of MADS-box gene family in Docynia delavayi (Franch.) Schneid. and its expression during different stages of seed germination, this study used seedlings at different stages of germination as materials and screened MADS-box transcription factors from the transcriptome database of D. delavayi using bioinformatics methods based on transcriptome sequencing. The physical and chemical properties, protein conservative motifs, phylogenetic evolution, and expression patterns of the MADS-box transcription factors were analyzed. Quantitative real-time PCR (qRT-PCR) was used to verify the expression of MADS-box gene family members during different stages of seed germination in D. delavayi. The results showed that 81 genes of MADS-box gene family were identified from the transcriptome data of D. delavayi, with the molecular weight distribution ranged of 6 211.34-173 512.77 Da and the theoretical isoelectric point ranged from 5.21 to 10.97. Phylogenetic analysis showed that the 81 genes could be divided into 15 subgroups, among which DdMADS27, DdMADS42, DdMADS45, DdMADS46, DdMADS53, DdMADS61, DdMADS76, DdMADS77 and DdMADS79 might be involved in the regulation of ovule development in D. delavayi. The combination of the transcriptome data and the qRT-PCR analysis results of D. delavayi seeds indicated that DdMADS25 and DdMADS42 might be involved in the regulation of seed development, and that DdMADS37 and DdMADS38 might have negative regulation effects on seed dormancy. Previous studies have reported that the MIKC* subgroup is mainly involved in regulating flower organ development. For the first time, we found that the transcription factors of the MIKC* subgroup exhibited a high expression level at the early stage of seed germination, so we speculated that the MIKC* subgroup played a regulatory role in the process of seed germination. To verify the accuracy of this speculation, we selected DdMADS60 and DdMADS75 from the MIKC* subgroup for qRT-PCR experiments, and the experimental results were consistent with the expression trend of transcriptome sequencing. This study provides a reference for further research on the biological function of D. delavayi MADS-box gene family from the perspective of molecular evolution.

Genome-wide identification of the banana GLR gene family and its expression analysis in response to low temperature and abscisic acid/methyl jasmonate / 生物工程学报

Glutamate receptor-like (GLR) is an important class of Ca2+ channel proteins, playing important roles in plant growth and development as well as in response to biotic and abiotic stresses. In this paper, we performed genome-wide identification of banana GLR gene family based on banana genomic data. Moreover, we analyzed the basic physicochemical properties, gene structure, conserved motifs, promoter cis-acting elements, evolutionary relationships, and used real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) to verify the expression patterns of some GLR family members under low temperature of 4 ℃ and different hormone treatments. The results showed that there were 19 MaGLR family members in Musa acuminata, 16 MbGLR family members in Musa balbisiana and 14 MiGLR family members in Musa itinerans. Most of the members were stable proteins and had signal peptides, all of them had 3-6 transmembrane structures. Prediction of subcellular localization indicated that all of them were localized on the plasma membrane and irregularly distributed on the chromosome. Phylogenetic analysis revealed that banana GLRs could be divided into 3 subclades. The results of promoter cis-acting elements and transcription factor binding site prediction showed that there were multiple hormone- and stress-related response elements and 18 TFBS in banana GLR. RT-qPCR analysis showed that MaGLR1.1 and MaGLR3.5 responded positively to low temperature stress and were significantly expressed in abscisic acid/methyl jasmonate treatments. In conclusion, the results of this study suggest that GLR, a highly conserved family of ion channels, may play an important role in the growth and development process and stress resistance of banana.

Identification and expression pattern analysis of RcACA gene family in castor under abiotic stresses / 生物工程学报

Auto-inhibited Ca2+-ATPase (ACA) is one of the Ca2+-ATPase subfamilies that plays an important role in maintaining Ca2+ concentration balance in plant cells. To explore the function and gene expression pattern of the RcACA gene family in castor, bioinformatics analysis was used to identify the members of the RcACA gene family in castor. The basic physical and chemical properties, subcellular location, protein secondary and tertiary structure, conserved domain, conserved motif, gene structure, chromosome location and collinear relationship, as well as the evolutionary characteristics and promoter cis-acting elements were predicted and analyzed. The expression pattern of the RcACA gene under abiotic stress was analyzed by expression (fragments per kilobase of exon model per million mapped fragments, FPKM) in castor transcriptome data. The results showed that 8 RcACA gene family members were identified in castor, acidic proteins located in the plasma membrane. In the secondary structure of all proteins, the α-helix and random coil is more; the RcACA genes were clustered into three categories, and the design of the genes in the same category was similar to the conserved motif. Both of them had four typical domains, RcACA3-RcACA8 had a Ca2+-ATPase N-terminal autoinhibitory domain. The RcACA gene is mostly located on the long arm of the chromosome and has 2 pairs of collinear relationships. There are more light response elements but fewer hormone-induced elements located upstream of the RcACA coding region. Interspecific clustering showed that the evolution of ACA genes among species was conservative. Tissue expression pattern analysis showed that RcACA genes showed apparent tissue expression specificity, and most of the genes showed the highest expression level in male flowers. Expression analysis under abiotic stress showed that RcACA2-RcACA8 were up-regulated under high salt and drought stress, and RcACA1 was up-regulated at 0-24 h under low-temperature stress, indicating that RcACA genes positively responded to abiotic stresses. The above results provide a theoretical basis for exploring the role of the RcACA gene in castor growth, development and stress response.

Genetic diversity analysis and fingerprints of Chrysanthemum×morifolium based on SSR molecular markers / 生物工程学报

The present study aims to explore the genetic diversity of germplasm resources of Chrysanthemum×morifolium (hereinafter, C.×morifolium) at the molecular level and to establish a fingerprint database of C.×morifolium varieties. We employed 12 pairs of primers with high levels of polymorphism, clear bands, and high degrees of reproducibility to analyze the SSR molecular markers and genetic diversity of 91 C.×morifolium materials and 14 chrysanthemum- related materials. With regard to constructing the fingerprints of the tested materials, we chose 9 pairs of core primers. The findings revealed that 12 primer pairs detected 104 alleles in 105 samples, ranging from 2 to 26. The average number of observed alleles (Na) per site was 9.25. The average number of effective alleles (Ne) per site was 2.745 6, with its range being 1.276 0 to 4.742 5. Shannon genetic diversity index (I) values ranged between 0.513 3 and 2.239 9 (M=1.209 0). Nei's gene diversity index (H) ranged between 0.216 3 and 0.789 1 (M=0.578 0). The observed heterozygosity (Ho) ranged between 0.223 3 and 0.895 2 (M=0.557 5). The expected heterozygosity (He) ranged between 0.217 4 and 0.793 3 (M=0.580 8). The polymorphism information content (PIC) ranged between 0.211 5 and 0.774 0 (M=0.532 9). The genetic similarity (GS) ranged between 0.228 5 and 1.000 0 (M=0.608 3). Cluster analysis revealed that when the genetic distance (GD) equals to 0.30, the tested materials can be classified into 2 groups. When the GD equals to 0.27, the first group can be divided into 6 subgroups; accordingly, 105 tested materials can be divided into 7 subgroups. The cophenetic correlation test was carried out based on the cluster analysis, and the corresponding results showed that the cluster map correlated with the genetic similarity coefficient (r=0.952 73). According to the results of Structure population analysis, we obtained the optimal population number, with the true number of populations (K) being 3 and the population being divided concerning Q≥0.5. Three subgroups, i.e., Q1, Q2 and Q3, included 34, 33 and 28 germplasms, respectively, and the remaining 10 germplasms were identified as the mixed population. During the experiment, 9 pairs of core primers were screened among the total of 12 for a complete differentiation regarding 105 tested materials, and the fingerprints of 91 C.×morifolium materials and 14 chrysanthemum-related materials were further constructed. Overall, there were significant genetic differences and rich genetic diversity among C.×morifolium materials, which would shed light on the garden application and variety selection fields of C.×morifolium. The fingerprint database of 105 C.×morifolium varieties and chrysanthemum-related species may provide technical support for future research regarding the identification and screening system of C.×morifolium varieties.

Gene cloning and enzymatic activity analysis of phenylalanine ammonia-lyase from Sinopodophyllum hexandrum (Royle) Ying / 生物工程学报

Phenylalanine ammonia-lyase (PAL) is the key entry enzyme of plant phenylpropanoid pathway. It plays an important role in the biosynthesis of podophyllotoxin, an anti-tumor lignan that is currently produced from its main natural source Sinopodophyllum hexandrum (Royle) Ying. In this study, we cloned the gene ShPAL encoding phenylalanine ammonia-lyase by RT-PCR from the root of S. hexandrum ecotype inhabited in the Aba' district, Sichuan, based on its public SRA transcriptome data-package. Bioinformatics analyses showed that the ShPAL-encoded protein is composed of 711 amino acids, contains the conserved domains of PAL, and has the signature motif within the active center of aromatic ammonia-lyases. Moreover, ShPAL protein was predicted to have a secondary structure mainly composed of α-helix and random coil, a typical 'seahorse' shape monomer tertiary structure, and a homologous tetramer three-dimensional structure by Swiss-Modelling. The phylogenetic lineage analysis indicated ShPAL was of the highest sequence identity and the shortest evolutionary distance with the PAL of Epimedium sagittatum from the same Berberidaceae family. Subcellular localization experiments showed that ShPAL protein was mainly distributed in the cytoplasm, despite of a minority on the endoplasmic reticulum membrane. Furthermore, ShPAL protein was recombinantly expressed in Escherichia coli and purified by histidine-tag affinity chromatography. Its enzymatic activity was determined up to 20.91 U/mg, with the optimum temperature of 41 ℃ and pH of 9.0. In contrast, the enzyme activity of its F130H mutant decreased by about 23.6%, yet with the same trends of change with temperature and pH, confirming that phenylalanine at this position does affect the substrate specificity of PAL. Both the wild type and the mutant have relatively poor thermostability, but good pH-stability. These results may help to further investigate the regulatory role of PAL in the process of podophyllotoxin biosynthesis and advance the heterologous synthesis of podophyllotoxin to protect the germplasm resource of S. hexandrum. They also demonstrate that ShPAL has a potential application in biochemical industry and biomedicine.

Functional analysis on sucrose transporters in sweet potato / 生物工程学报

Sweet potato is an important food crop that can also be used as an industrial raw material. Sucrose is the main form of long-distance carbohydrate transport in plants, and sucrose transporter (SUT) regulates the transmembrane transport and distribution of sucrose during plant growth and metabolism. Moreover, SUT plays a key role in phloem mediated source-to-sink sucrose transport and physiological activities, supplying sucrose for the sink tissues. In this study, the full-length cDNA sequences of IbSUT62788 and IbSUT81616 were obtained by rapid amplification of cDNA ends (RACE) cloning according to the transcripts of the two SUT coding genes which were differentially expressed in sweet potato storage roots with different starch properties. Phylogenetic analysis was performed to clarify the classification of IbSUT62788 and IbSUT81616. The subcellular localization of IbSUT62788 and IbSUT81616 was determined by transient expression in Nicotiana benthamiana. The function of IbSUT62788 and IbSUT81616 in sucrose and hexose absorption and transport was identified using yeast functional complementarity system. The expression pattern of IbSUT62788 and IbSUT81616 in sweet potato organs were analyzed by real-time fluorescence quantitative PCR (RT-qPCR). Arabidopsis plants heterologous expressing IbSUT62788 and IbSUT81616 genes were obtained using floral dip method. The differences in starch and sugar contents between transgenic and wild-type Arabidopsis were compared. The results showed IbSUT62788 and IbSUT81616 encoded SUT proteins with a length of 505 and 521 amino acids, respectively, and both proteins belonged to the SUT1 subfamily. IbSUT62788 and IbSUT81616 were located in the cell membrane and were able to transport sucrose, glucose and fructose in the yeast system. In addition, IbSUT62788 was also able to transport mannose. The expression of IbSUT62788 was higher in leaves, lateral branches and main stems, and the expression of IbSUT81616 was higher in lateral branches, stems and storage roots. After IbSUT62788 and IbSUT81616 were heterologously expressed in Arabidopsis, the plants grew normally, but the biomass increased. The heterologous expression of IbSUT62788 increased the soluble sugar content, leaf size and 1 000-seed weight of Arabidopsis plants. Heterologous expression of IbSUT81616 increased starch accumulation in leaves and root tips and 1 000-seed weight of seeds, but decreased soluble sugar content. The results obtained in this study showed that IbSUT62788 and IbSUT81616 might be important genes regulating sucrose and sugar content traits in sweet potato. They might carry out physiological functions on cell membrane, such as transmembrane transport of sucrose, sucrose into and out of sink tissue, as well as transport and unloading of sucrose into phloem. The changes in traits result from their heterologous expression in Arabidopsis indicates their potential in improving the yield of other plants or crops. The results obtained in this study provide important information for revealing the functions of IbSUT62788 and IbSUT81616 in starch and glucose metabolism and formation mechanism of important quality traits in sweet potato.

Identification, expression and DNA variation analysis of high affinity nitrate transporter NRT2/3 gene family in Sorghum bicolor / 生物工程学报

Nitrate is the main form of inorganic nitrogen that crop absorbs, and nitrate transporter 2 (NRT2) is a high affinity transporter using nitrate as a specific substrate. When the available nitrate is limited, the high affinity transport systems are activated and play an important role in the process of nitrate absorption and transport. Most NRT2 cannot transport nitrates alone and require the assistance of a helper protein belonging to nitrate assimilation related family (NAR2) to complete the absorption or transport of nitrates. Crop nitrogen utilization efficiency is affected by environmental conditions, and there are differences between varieties, so it is of great significance to develop varieties with high nitrogen utilization efficiency. Sorghum bicolor has high stress tolerance and is more efficient in soil nitrogen uptake and utilization. The S. bicolor genome database was scanned to systematically analyze the gene structure, chromosomal localization, physicochemical properties, secondary structure and transmembrane domain, signal peptide and subcellular localization, promoter region cis-acting elements, phylogenetic evolution, single nucleotide polymorphism (SNP) recognition and annotation, and selection pressure of the gene family members. Through bioinformatics analysis, 5 NRT2 gene members (designated as SbNRT2-1a, SbNRT2-1b, SbNRT2-2, SbNRT2-3, and SbNRT2-4) and 2 NAR2 gene members (designated as SbNRT3-1 and SbNRT3-2) were identified, the number of which was less than that of foxtail millet. SbNRT2/3 were distributed on 3 chromosomes, and could be divided into four subfamilies. The genetic structure of the same subfamilies was highly similar. The average value of SbNRT2/3 hydrophilicity was positive, indicating that they were all hydrophobic proteins, whereas α-helix and random coil accounted for more than 70% of the total secondary structure. Subcellular localization occurred on plasma membrane, where SbNRT2 proteins did not contain signal peptides, but SbNRT3 proteins contained signal peptides. Further analysis revealed that the number of transmembrane domains of the SbNRT2s family members was greater than 10, while that of the SbNRT3s were 2. There was a close collinearity between NRT2/3s of S. bicolor and Zea mays. Protein domains analysis showed the presence of MFS_1 and NAR2 protein domains, which supported executing high affinity nitrate transport. Phylogenetic tree analysis showed that SbNRT2/3 were more closely related to those of Z. mays and Setaria italic. Analysis of gene promoter cis-acting elements indicated that the promoter region of SbNRT2/3 had several plant hormones and stress response elements, which might respond to growth and environmental cues. Gene expression heat map showed that SbNRT2-3 and SbNRT3-1 were induced by nitrate in the root and stem, respectively, and SbNRT2-4 and SbNRT2-3 were induced by low nitrogen in the root and stem. Non-synonymous SNP variants were found in SbNRT2-4 and SbNRT2-1a. Selection pressure analysis showed that the SbNRT2/3 were subject to purification and selection during evolution. The expression of SbNRT2/3 gene and the effect of aphid infection were consistent with the expression analysis results of genes in different tissues, and SbNRT2-1b and SbNRT3-1 were significantly expressed in the roots of aphid lines 5-27sug, and the expression levels of SbNRT2-3, SbNRT2-4 and SbNRT3-2 were significantly reduced in sorghum aphid infested leaves. Overall, genome-wide identification, expression and DNA variation analysis of NRT2/3 gene family of Sorghum bicolor provided a basis for elucidating the high efficiency of sorghum in nitrogen utilization.

The role of plant WRKY transcription factors against salt stress: a review / 生物工程学报

High salt content in soils severely hampers plant growth and crop yields. Many transcription factors in plants play important roles in responding to various stresses, but their molecular mechanisms remain unclear. WRKY transcription factors are one of the largest families of transcription factors in higher plants that are involved in and influence many aspects of plant growth and development. They play important roles in responding to salt stress. The regulation of gene expression by WRKY proteins is mainly achieved by binding to the DNA's specific cis-regulatory elements, the W-box elements (TTGACC). In recent years, there have been many studies revealing the roles and mechanisms of WRKY family members, from model plant Arabidopsis to agricultural crops. This paper reviews the latest research progress on WRKY transcription factors in response to salt stress and discusses the current challenges and future perspectives of WRKY transcription factor research.